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Abstract: Fusion reactions with a weakly bound projectile are studied using the double-folding model along with 

a repulsive interaction modifying term. Using this modified potential, including nuclear matter incompressibility 

effects, the fusion reaction cross sections and suppression parameters are calculated for ^Be + 

and Al reactions. The results show that applying these effects at energies near the Coulomb barrier improves the 

agreement between the calculated and experimental cross sections, and modifies the mean values of the suppression 

parameter. 
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1 Introduction 

In recent years, there have been many theoretical and 
experimental studies of heavy-ion reactions with light 
projectiles which have low binding energy [1-5]. Such 
weakly bound nuclei as and ®i?e are of great impor¬ 
tance in forming super-heavy nuclei and understanding 
the nuclear structures. As a result of the low binding en¬ 
ergy of such nuclei compared with tightly bound nuclei, 
it is interesting to follow the reaction through various 
channels including complete fusion (CF) and incomplete 
fusion (IGF). It has recently been reported that the com¬ 
plete fusion cross sections of such reactions in theoretical 
calculations where breakup is not included are about 25- 
30 percent greater than the experimental cross sections 
[6]. The absolute value of this difference is known as 
the suppression parameter. This parameter in the com¬ 
plete fusion cross section seems to originate from parallel 
channels of the reaction during the CF direct process, 
which competes with the main channel of the expected 
reaction. Studies in recent years have shown that the 
input channels of a reaction, such as total interaction 
potential, have a significant impact on changing this pa¬ 
rameter. This study aims to investigate the change of 
rate of this parameter as a result of the application of 
modifying effects to the double-folding potential in ^Be 
reactions with ’^^Si and as target nu¬ 

clei, the experimental data of which have recently been 
presented [7-10]. These reactions are studied using the 
modified double-folding potential and coupled-channel 
method to calculate the complete cross sections. We 


also want to compare light target nuclei such as 
and with ‘^^^Bi, and explore whether or not 

complete fusion is also affected by the breakup process. 
This is an extension of our previous work [II] in which 
we conclude the trend of complete fusion suppression in 
reactions with light and heavy targets. 

This paper is organized as follows: a description of the 
theory employed to calculate the nuclear potentials by 
using the double folding (DF) model and simulation of 
incompressibility effects is given in Sec 2. In Sec. 3, we 
present and discuss the effects of these corrections on 
the calculations of the fusion cross-sections. Section 4 
gives some concluding remarks. 

2 Model for ®Be-induced interaction potential 
2.1 Basic interaction potential 

One of the affecting factors in heavy-ion fusion reactions 
is the total interaction potential. Since the Coulomb and 
angular parts of the potential can be calculated with a 
high level of accuracy, there remains only the nuclear 
potential to be calculated in each interaction. There are 
several methods for calculating the interaction poten¬ 
tial. Among these models, the proximity type [12-16] 
and Sao-Paulo [17] potentials have been mostly used for 
various reactions. The phenomenological base of these 
models motivated us to calculate the real part of the to¬ 
tal interaction potential with the double-folding model 
[18]. The following integral form is used to calculate the 
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nuclear potential: 


V{r) 


drp 


drTPp{rp)pT{rT)vi2{s) 


( 1 ) 


where s = r+rp—r^. In this relation Pt{yt): pp(vp) and 
Vi 2 [s) are nucleon densities of the target and projectile 
and the central part of the NN interaction between two 
nuclei respectively. 

To reproduce the G-matrix elements of the NN poten¬ 
tials in an oscillator basis the M3Y-Paris formalism is 
a popular choice [18]. This density-independent M3Y 
interaction has been used with some success in folding 
model calculations of the heavy-ion potential, and its 
explicit form follows below. The direct part is 

(^) = - 25372 Eyi! 2 (MeC) 

( 2 ) 

and the zero-range pseudo-potential can be used reliably 
to represent the knock-on exchange part. The Paris pa- 
rameterizations of zero-range approximation for the ex¬ 
change part of the interaction can be expressed as [15] : 

= -590(1 - 0.002 ^)(Mel/) (3) 


where ^ is the bombarding energy per projectile nucleon. 
The experimental density distributions are used both for 
target and projectile. The density distribution function 
of the projectile nucleus ®Be is assumed to be of the the 
Harmonic Oscillator form: 


pp{r)=po 


1 + Q!p(- 

dp 


exp 



( 4 ) 


where the coefficients ap and Op are 0.611 and 1.791 
respectively [15]. This form of density profile offers an 
effective spherical shape for the ground state of ®Be, so 
it is possible to neglect the deformation effects of the 
nucleus and the rotational band of the ground state of 
®Be to calculate the nuclear potential. For parameteriza- 
tions of target densities, we have used the two-parameter 
Fermi model (2PF) [19] . 


P2PF{f) 


Po 


l-bexp 1^- 


( 5 ) 


The radius i? 2 pp and the diffuseness a 2 PF in the above 
equation are listed in Table 1. 


Table 1. Experimental density distribution pa¬ 
rameters for target nuclei |19[ I20| . 


Nucleus 

Ropf (fm) 

a2PF(fm) 


6.75 

0.468 

206 pb 

6.631 

0.505 


3.017 

0.52 


3.07 

0.519 


The six-dimensional integral, Eq. 1, can be calcu¬ 
lated using Fourier transformation in momentum space. 
Thereby one can reduce this integral to a product of three 
one-dimensional integrals using 



drpdrTPp{rp)pT{rT) 


VM3Y{s)e^ 


= J y drpdrppp(rp)/9p(rT) J UM3v(s)e' 7K.(s+rj, 

= y/9p(rp)e“‘‘"'’^drp J pTirT)e~"'"’"^drT J ■Umsv 

= pp(k)pp(k)p(k) (6) 

In these calculations the Fourier transformation is 
given by: 

poo 

f{k)=47r f{r)jo{kr)r^dr (7) 

and its inverse transform is: 

poo 

/(r) = l/(27r®) / f{k)jo{kr)k'^dk (8) 

Jo 

In the calculation of total potentials the Coulomb po¬ 
tential is 


K = (2’pZpeV2i?)(3 - rVi?®) r < i? 

= iZpZTey2R)r>R (9) 

where R = 1.2^4^^^ and Ap is the mass number of the 
target nucleus. 

The process of density distributions overlapping of 
interaction nuclei implies that the variation of total 
density of colliding nuclei can be increased up to the 
value of the saturation density of nuclear matter i.e. 
Po — 0.17/m“®. This non-physical effect can be modified 
by including the effect of equation of state of nuclear 
matter in calculation of the nuclear potential. This is 
explained in the following. 

2.2 Modified interaction potentials 

We expect that the interaction potential between tar¬ 
get and projectile nuclei will contain an additional re¬ 
pulsive core in the interior region based on the Pauli 
Exclusion Principle, which prevents the complete over¬ 
lapping of the wave functions of two fermionic systems 
[21]. The DF model is one of the theoretical models 
which uses the sudden approximation to analyze the fu¬ 
sion process. Therefore, based on this model, the total 
density in the overlapping region can be increased to val¬ 
ues up to p ~ 2po. Therefore this increasing energy in 
the total potential at this region can be estimated by 

AU^2Ap[e{2po,S)-e{po,S)] (10) 
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where Ap is the mass number of the projectile nucleus, 
5 is the relative neutron excess and e(po,S) is the energy 
per nucleon as evaluated by the Thomas-Fermi equation 
of state for cold nuclear matter [21]. By employing a 
repulsive core in the N-N interaction, Vrep(r) = 

, one can simulate this increasing energy in the DF cal¬ 
culation. 


Vl2(r)=Vatr(r)+Vrep(r) (11) 


where Vatr(r) is the attractive part of the NN interaction 
and Vrep(r) is the repulsive core which appears from the 
Pauli exclusion principle. We have used the DF method, 
Eq.(l-8), to calculate the attractive and repulsive terms 
of the nuclear part of the total potential. In calculating 
the repulsive part of the nuclear potential, it is assumed 
that the diffuseness parameter for the density distribu¬ 
tions of target and projectile are equal to a^ep- The free 
parameters a^ep and z/^ep are adjusted so as to meet the 
following condition. 



r(fm) 



r(fm) 



r(fm) 


= ( 12 ) 


where Kep(O) is equal to the repulsive part of the total 
potential evaluated at r = 0 by the DF integral, and K 
is the incompressibility of cold nuclear matter. 

3 Calculations and resnlts 

In this study the total interaction potentials have 
been calculated in momentum space using Eq. (1-8) 
for ®Be-|-^°®Bi,^°®Pb,^®Si,^’^Al fusion reactions. Harmonic 
oscillator (HO) and two-parameter Fermi (2PF) density 
distributions have been used in the calculations for pro¬ 
jectile and target nucleus distribution functions. These 
interactions are named M3Y and depicted in Fig. 1 as 
dashed lines. Applying the modifications related to the 
nuclear matter saturation and adding the repulsion core 
term to the M3Y potential, the interactions are named 
MIP and shown by solid lines in Fig. 1 . 


Fig. 1. Radial shape of the folded potentials 
for the ®Be-h^°®Bi,2°®Pb,29Si,®'^Al fusion reac¬ 
tions calculated using the M3Y-Paris interaction 
(dashed lines) and MIP results (solid lines). 


As can be seen in Fig. 1 the inner regions of the M3Y po¬ 
tential include a very deep valley, meaning a very strong 
absorption in that region which is physically impossible. 
With the applied modifications, a valley of finite depth 
emerges, which seems to be necessary in the known M3Y 
potential. The obtained results for the interaction po¬ 
tentials were used to calculate the fusion cross sections 
using the CCFULL coupled-channels code [22]. The cal¬ 
culated cross sections based on M3Y and MIP interac¬ 
tions are compared to experimental data for the fusion 
of ®Be-h^°®Bi,208pb,29Si,27Al in Fig. 2. 
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Fig. 2. The fusion cross sections for the 
®Be+^°®Bi,^°®Pb,^®Si,^^Al systems are compared 
to the couple-channel calculations using M3Y in¬ 
teraction (dashed line) and MIP interaction (solid 
line). The experimental data also have shown by 
solid circles [7-10]. 


In Fig. 2 the dashed lines indicate the M3Y based 
cross sections and the solid lines indicate the MIP based 
cross sections. The calculated cross sections using MIP 
interactions, which are based on nuclear matter incom¬ 
pressibility, show a better agreement with experiment 
than calculations based on the M3Y interaction. This 
was achieved by adding the repulsive part of the NN 
interaction that is used in the calculation of DF nuclear 
potentials. From these figures the existence of adding a 
repulsive core to the calculation of interaction potential 
seems to be helpful in describing weakly bound fusion 
systems. 
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Fig. 3. Suppression parameter for the 

®Be-t^°®Bi,2°®Pb,2^Si,2'^Al fusion systems. The 
solid circles indicate the M3Y based calculations 
and the open circles indicate the MIP based cal¬ 
culations. The mean values for each reaction are 
indicated in the figure. 


A wide variety of studies on weakly bound nuclei re¬ 
actions in recent years have introduced a modifying term 
[23,24]. These investigations suggest the necessity of in¬ 
troducing the suppression parameter (SP) in the theo¬ 
retical calculations. This parameter is dependent on the 
potential model used in the interaction and is defined as 
follows 

5'F’= (1 — crBjjp./CTrfieor.) X 100 (13) 
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where (Texp. and cfTheor. are fusion cross sections ob¬ 
tained from the experimental measurements and theoret¬ 
ical calculations, respectively. Therefore, the repulsion 
term effect resulting from overlapping nuclear densities 
of the two interactive nuclei on the suppression parame¬ 
ter mean value is calculated using Eq. (13) for the MIP 
and M3Y models and displayed in Fig. 3. in terms of 
Ec.ruJ^B- 

It can be observed from this figure that the nuclear 
matter saturation effects can modify the suppression pa¬ 
rameter values. The changes are listed in Table 2. 

Table 2. Suppression parameter values for fusion 
reactions before and after the compressibility of 
nuclear matter is included. 


Nucleus 

SP(M3Y) 

SP(MIP) 

209 Bi 

18.27 

6.14 

208 pjy 

46.00 

44.98 


22.46 

20.20 


14.31 

11.28 


The differences which we see in Table 2 for the SP(M3Y) 
in the ®Be-|-^°®Bi,^°®Pb, systems is 28 percent. Fusion 
reactions with ®Be,®’^Li beams incident on ^°®Bi and 
^°®Pb targets have investigated in [25-31]. Above-barrier 
CF cross sections for the reaction of ®Li with ^°®Bi and 
^°®Pb are in agreement with each other and both show 
about 34 percent for the SP value [28,29,32]. How¬ 
ever, the situation is not so clear for the ®Be-induced 
reactions, where the difference between measured com¬ 
plete fusion cross sections for the ®Be-|-^°®Bi system and 
the neighboring ®Be-|-^°®Pb system is about 25 percent 
[33-36]. Furthermore, discrepancies exist even between 
the different sets of complete fusion measurements for 
®Be-|-^°®Bi, which were carried out in different labora¬ 
tories [27,37]. The large difference in above-barrier CF 
cross sections between reactions with ^“®Bi and ^°®Pb 
targets is discussed in [38], which attributes the current 
problem to normalization for one or more measurements. 

4 CONCLUSIONS 

In summary, we have modified explicitly some gener¬ 
alized and realistic calculation processes for the double 
folding model including nuclear matter incompressibil¬ 
ity. The modified program is applied to calculate the 
real part of the interaction potential for several weakly 
bound fusion systems. This study also examines the 
effects of modeling the repulsion term in the nuclear po¬ 
tential on the fusion cross sections and the effect of the 
suppression parameter of nuclear fusion reactions with a 
weakly bound projectile. The studies have been carried 


out in this case on ^Be + and 

reactions. The results show that applying the nuclear 
matter incompressibility effects to the nuclear potential 
calculation limits the nuclear potential depth in the in¬ 
ner regions of the ion-ion distance. As shown in Fig. 2, 
this improves the agreement between the calculated and 
experimental cross sections. The applied modifications 
also modify the mean value of the suppression parame¬ 
ter, which is displayed in Fig. 3 for ^Be + '^°^Bi, 

‘^^Si and reactions. 
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